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CHAPTER 1 


INTRODUCTION 


Let xo xe ae ee b a i i } 

( 1 Y4)» ( 5 Yen) . Gy) e a sample from a bivariate 
normal population with correlation po. If either (x) 5%5, etenses x) or 

(Yp Yo» Noon J.) is a series of independent observations, the sample product- 
moment correlation r has the familiar null distribution (p=0), (see 
Keeping [7]), but, if both series are serially correlated, very little is 
known about the distribution of r, Orcutt and James [12] have made the 
following remarks relevant to this problem: 

The testing of the significance of a correlation involves 

a comparison with what would have been obtained between 

non-related series thought to be analogous to the observed 

SCGMCS uae. 6 = 

There is, however, one obvious point at which the analogy 

underlying such tests of significance may break down when 

one is concerned with economic time series, namely, if 

the consecutive terms are really correlated. In economic 

time series, in metereological time or spatial series, or, 

for that matter, in biological time series, auto-correlation 

usually exists. ..... 

In economics, most of the material that we wish to investigate 
for relationships exhibits auto-correlation, and there is 

a real need for a test of significance for correlations 

which is based on a more realistic sampling model. 

The literature pertaining to this subject is sparse, For a concise 
review of some of the earlier studies of the distribution of r , we quote 
McGregor's paper [9]: 

When both series are of the stationary Markov type with 


auto-correlations p, and Pp » Bartlett (1935), [2], 


has shown that the approximate variance of r is given 
by 


“touboxg ai genes 


+ 


» Sehenoa ® ren 


PO 


Deets 
var(r) ~ AOC Sas Bt : 
ee 


Orcutt and James (1948), [12], concluded, from the results 
of a sampling experiment, that there is some evidence that 
the variance of the correlation between two time-series is 
nearly independent of the true auto-correlations, In another 
sampling experiment, Quenouille (1949), [15], worked out 
various partial cross-correlations between two time-series 


with P, = 0.6 and , = O.4. In particular, the observed 


variance of r(x,y,/%,¥,); the partial cross-correlation 


between x, and Y, with the effects ey and Yoo 

removed, was considerably less than that of r, Hannan 

(1955), [6], has shown that r(x,¥51%,¥,) provides a more 
(= 


efficient test of 090 =0 than does fr, 


Using an approach due to Daniels [3], [4] based on the method 
of steepest descents, McGregor [9] obtained the approximate null distribu- 
tion of r, when the series (x,] and Wei are each of the stationary 


and 90., 


linear Markov type with know means and auto-correlations Py 2 


respectively, 


In this thesis we shall consider the more general case where the 
means are fitted, and derive the approximate null distribution of r _ for 
the two stationary, linear Markov processes (x, } and fy,) with serial 


correlations Py and Pos respectively. 


As conjectured by McGregor, the analysis of [ 9] can be extended 
to the fitted means case, In Chapter II we discuss McGregor's method for 
the known means case in detail in order to keep clarity and unity throughout 
the development of the more general case where the means are fitted, The 
approximate density function (p*(r)) of the product-moment correlation r 
for the fitted means case (as for the known means case [9 ]) was fo! to 


depend on the auto-correlations Py and Po only through their product, 
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[see (3.33) (also ( 2.61))]. Upon renormalization of p*(r), (3.46), we 
observe that it is of the same form as p(r), the approximate density 
function of r for the known means case, with N replaced by (M-1). The 
main problem in this development was the approximation of the determinant 
|B | encountered in the joint-moment generating function of the quadratic 


forms used in defining r, (see Appendix IIIB). 
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CHAPTER II 
THE APPROXIMATE DISTRIBUTION OF THE CORRELATION 


BETWEEN TWO STATIONARY, LINEAR MARKOV SERIES WITH KNOWN MEANS 


In this chapter we discuss in detail McGregor's [9] derivation 
of the approximate null distribution of r when the series {x, } and (ya! 
are each of the stationary, linear Markov type with known means and with 
and 0 


autocorrelations 0p respectively. 


il en 


Consider the two stationary, linear Markov processes 


x, = 0, X + €, for all i 
and 


Feel + 1; for all...3.; 


where the e's and 1's are independent N(0,1) variables, The joint dis~ 


tribution of Xp Xo» seen KX 


2 x 4 5 1 i « pod 
dF, (x).--->%,) = (2n) iz | 2 exp E cope x féx..-dy,, 


(axn) ; 
where! 9x )— (X)5%52-+-9%))! and Le / is the covariance matrix of x. 


For the stationary, linear Markov process it is easily shown that 


A 
var (x, ) =1- 0° 
and lc | 
ral 
cov(x, | %;) = ane (for any fixed k) . 


Hence, the covariance matrix for x is 
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Also, as in Patton's thesis [13], it can be shown that 
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and hence, 


Q els 2 
x i. x= sr} Py 20, (x, 


Similarly, the joint distribution of VprYorerery 


dF (yj >+++2¥,) = (2m) 


n-2 
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2 
1-9, 
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where y = (¥) Yor seo %n) and aaa is the covariance matrix of y 


which is of the same form as (2.1) except that P1 is replaced by Po: Also, 
bea and z* are given by (2.2) and (2.3), respectively, with P, 


= 
replaced by Po Hence, 
1 =] 2 
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Now the joint distribution of Kp Kore skye VyrVore eV, is 
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We seek an approximate distribution of the product-moment correla- 


tion 


where 
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The joint moment-generating function for the distribution of C, D and E is 
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The exponent of the integrand is 


TC+SD+UE - ax Piha ihe x| 


ane ax fo z+ yy to y - 2lC = esp - ave | 
* eye 


mee Dae 2 2 p02 ne 2 2 
= 3 [xi + (1407) (x5 NAA Yi) See en) eae yaa (1405) (yo+-+-+¥, 1) tye 


- 20, (xX +--+ 1X) - 207 ete t95 yd 


2 2 2 2 
- aT(xi+...+x) - a8 (yj+..-+¥_) - 2U(x,y,+-- +,7,) | 


i 
i] 
[) 
_- 
1 
‘|x 
Rw 
> 
wo 
ALE 
ww 


ii 
i 
DY) og 
* 
x 
[> 
*% 
* 


TS Onset erh see TOi, wohtogid gr Javeca: 
‘ 7m iN (oO. ae ai hy 


H 


‘ 


x t ql 
where x* = oy == (Xp 520s oVyreee sy ) and A is the en x en 


partitioned matrix, 
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and Qs the scalar matrix with diagonal elements (-U). 


Thus (2.5) is 


2 ay 2 oy 
(Ize;)*(i-05)F pe 00 
M2), 0) = Hse i exp[-3 x*'A x*] 
(2) “ee9 ios 
xdx,...dx dy,...dy, 5 
Since A is a symmetric,nonsingular 2n x 2n matrix and x*' = (Xp oee69%o 


ES A there exists an orthogonal transformation of the variables 


ules nig ‘ 


MO us 
mh Ne Be) Hb 


be 


cep 


ie 


a i 


} ab ‘er pe } 


| for : 


x acre XY, AOC; which diagonalizes A, Making this transformation and 
s ie é 
performing the integration, we get 
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M(T,S,U) = ———"—_,-—* ,, (2.6) 
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where |A| is the Jacobian of the transformation. 


In Appendix III, (1II.23), the determinant of A is evaluated 
approximately giving 
2 2 
[a5(By4745) = 25(By9781) (Bp) -8)! ie 
n 2 . ; 
a,(1-a)+a,,)(1+a,+a,)(1-a,) 
where terms which will ultimately become exponentially small have been 


omitted, and 
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From (2.6) and (2.7) we obtain 
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Following the development of Appendix I, we make the substitution 


me 


v(r +r, z : 
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’ (2.11) 


Shand = oe gal Soe = . By the relations (2.8) and (2.11) a, and a 


Eee 2 1 2 
are implicit functions of T, S and vY (a, 2a,(T,3,v), a, = a,(T,8,v)) and 
hence 
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where the error is relatively o(+) : 


Then using the inversion formula (1.3), we obtain the following joint dens- 


tity ion 35 and r_! 
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where terms which are relatively o(=) have been ignored. 
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(Note that there is a distinction between a, = a(T,S,v) and 


a(T,S), similarly between a, = a, (T,S,v) and 


al yao 
arleee = a,(T,8), but throughout the discussion it will be 
apparent which of these functions is used.) 

Now putting v = 0, so that a, and ay are functions of T and 


S only, and eliminating a, from equations (2.8), we may write 
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Similarly, differentiating a, with respect to S (holding T constant) 
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From (2.15) and (2.16) we solve for T and § simultaneously as follows: 
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Finally from (2.15) and (2.16), solving for $, we find 
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Then 


A A Ane (2 ee 
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With this notation and using (2.18), we have 


2 A 
1+p,-er i 1 
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- 84, r® i 51 Pa (1405) + 2(1+48 o)e5e 1 pp 140, ) 
a2 2 ; ae, 2.2 2; 
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Now substituting (2.22), (2.23) and (2.24) in the above expression , we get 
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Substituting the above expressions in a(B54 a) - a5(Bi5 a, (Bo, a,)> we 
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Differentiating (2.35) partially with respect to w and holding T and 
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approximate distribution reduces to the familiar null distribution of the 
product-moment correlation with known means, which is 
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CHAPTER III 
THE APPROXIMATE DISTRIBUTION OF THE 


CORRELATION BETWEEN TWO STATIONARY, LINEAR MARKOV 


SERIES WITH FITTED MEANS 


In this chapter we derive the approximate null-distribution of the 
sample product-moment correlation, r, when the series {x,} and {y,) 
are each of the linear, stationary Markov type where the means are fitted 
and the auto-correlations are Py and Pos respectively. In this deriv- 
ation, which follows the general stages of the procedure in Chapter II, we 
encounter the matrix B (3.6) whose determinant, |B| , is evaluated 


approximately in Appendix III. 
Consider the two stationary, linear Markov series 
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where the error is relatively o(=) : 


Then using the inversion formula (1.3), the joint density of ty and t, is 
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where the terms which are relatively 0(=) have been omitted, 
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To evaluate the integral 


5 3/2 fa) | 
PP a Bo -5) (1- “a, a) (l+a 1+ 7 “(lea >) =|" ioe 
J a x - -a \ 
as J ae a,)-a5(B1572,)(Bo1-24)1 dis 
(3.13) 
we use the bivariate, saddle-point approximation, 
ee Pa (W(2, 52, ne @ nen )dzdz, 
* k™ i 
9'(241%,) (W(2,,2,))* * 
, (3.14) 


pF ie ee ey ale a 
2xk™{Y, (2) 02, Wool 2,025)" [¥10(2, 925) 173? 


A A A A A 1), \ 
where ¥1,(2) 925); Voon( 24225) Wyo(2, 225), 2, and Z, are defined in (2.14). 


In (3.13) T and S correspond to Zz, and Z, in (3.14), respectively, 


Slag to W(24525), 3(n=5) to k and 


pS) Ec IA ay to o”(z 22 ) ; 
(2(8yy78)-85(81078,)(Boq7@,)]1 1% ee 


(1-a,42,)2(1+a,4a,)°/*(1-a,) 3 2 


v=0 


The paths of integration in the T and S_ planes are taken as the lines 


of steepest descent of a 3 a(T,S) such that 


al veo 
0a,(T,8) 1 da,(T,S) 

of 38 
Now comparing (3.13) and (3,14) with (2,13) and (2,14), it can be seen that 


@'(T,8) = (1+a,+a,) 9(7,8) 


and 


Sy rk bend? 36 sxe 
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= #(n - 5), while k = S(t 
(note: by replacing n by (n-1) in k we obtain 


k*) 


Otherwise there are no differences between (3.13) and (2.13), and hence to 
* A 
obtain (*,8) ; 


* A A A 
@ (2,8) = (144,+4,) 9(f,8) 
A AVN A 
where a, = a, (T,S) and a 


a(T 59), we can simply use the expressions 
(2.34) and (2.30) derived in Chapter II 
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= Sec. = @ (ys) 

[a (B,174,)-85(8, 74, ) (8-8 a_)] as 
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From (2.40), we have 
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Finally, since * = 3(n-5) and V(2) 525) = a(T,S), the bivariate saddle- 


point approximation of (3,13) is 
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Thus the approximate joint marginal density function of aT and ty 
(3.12) becomes 
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A 1 DU SAO UNOS 
n*(t) 575) 3(n-1 2 
(1-9, )(1-0,)(9,°,) (r,+r,) 
2 3/2 mS N74 ANE A Sa-2 
_ oeeittad eatenta) » (t-ea)e(tte2)) 25 


n(n-5)(1-a,,0,0,,) 
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where the relative error of the approximation is 0(4) 


The product-moment correlation r is related to T and r,s 


r The approximate density function of r can be obtained from 


= Sa Saye 
(3.16) as a marginal density. First let 


C = tT, and €& = A ‘ 
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Writing the expressions involving ty and ty in (3.16) in terms of ¢ 


and €, we have 
(+, (2+0,)"44 (1-0 ,)°) = Ve [Ve (1-05) + Fp (ee) ] 


[x (1405)+r( 1404) ] = Ve [ vé (1405) + € (1405)| ; 


[x (1-05)+r,(1-07)] =e [ve (1-05) + i: (1-04) ]] 


and 


by 


ahs 


(Bree Sieg Tera! eye ge! 


em em a 


a 
pean 
eas in 
ere 


= 6% «= 


'p) nD °n9Y 5 ‘ ) 1 
{x5 (1-05)"4n5 AGS p' =) 42 x al (e,-e,) “+ (1-p,0,)7] i 
i Pond : ag } 


sac i “a 2 be Ne ) +2[(0 5-05)” *+(1- “P1P5) ae 


Thus, substituting the above expressions in (3,16), the approximate joint 


density of (€ and € is 
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(3.17) 
where a, is now a function of € and €., 
Then the approximate marginal density function of @ is 
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Following Daniels [], the saddle-point approximation of (3.18) is 
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Using the relations (2.59) and (2.60), we may write (3.32) as 
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Since the approximate distribution depends only on the product of the serial 
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where the error of the approximation is relatively o(+). We note that the 
right member of (3.33) is an even function of r and hence to the order 


considered, the mean and all of the odd moments are zero, 


From the approximate density function p’(r) (3.33), we may find 


the approximation of the variance of r. From (2.62) we can write 
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Denoting the integral of (3.34) by I and using (2.64), we may write, 
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Therefore, the variance of r is approximately 
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where the error is relatively o(+). Thus, using the density function 
p*(r), which is derived in this chapter, we have verified Bartlett's [2] 


approximation of the variance of r. 


To renormalize the approximate density function p*(r), we 
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Thus, let us adjust the exponents of (4.37) so that 
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As in Chapter II, we expand the following functions of : as Maclaurin 


series, in order to obtain the most convenient 5 3 


for g,(r°) = (Giepethere le Lhuddne (2.69), ‘we have 

g,(r*) = (14a)"/? cok mor’), (3-39) 
for g,(r°) s [t( 140)? wor? + (140) 77? , using (2.70), we have 

gn(r*) = [2(14a)]"9/? hax2] PCa + o(r') (3.40) 
and for g(r") = [t (140) odour? 2 + (1-0)]° , using (2.71), we have 

g(r ) = 2 fue] Oe + o(r") i (3.41) 


Since the expected value of n° is o(n”*), the error involved in substi- 
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For the renormalization of p*(r) we must evaluate the integral 
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sample product-moment correlation r depends only on the product of 
serial correlations, Py and Po» and that for a = PP, = 0, p (xr), 
(5.46), reduces to the familiar null distribution of the product-moment 


correlation with fitted means, (see Keeping (11.12.3) in [ 7] 
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We also note that by replacing N by (M-1) in p(r), (2.76), the approx- 
imate density function of r for the known-means case, we obtain the approx- 


imate density function of r for the fitted-means case, namely, P(r), 
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For a graphical representation of p*(r) (3.46), see Figure 1 


below. For the calculations Peters' [14] and Davis'[5] tables were used, 
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Figure 1, 
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p(x) (approx. ) 


The approximate density function of the product-moment correla- 
tion r between two stationary linear Markov series with 
fitted means for parameter values a = -0.5, 0.0, 0.5 amd for 
sample size n = 30. 
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APPENDIX I 


A METHOD OF INVERSION 


It will be convenient to consider here the method of inversion, 
developed by McGregor [9], for determining the joint distribution of two 


ratios with a common numerator. 


Let £(C,D,E) be the probability density function of the joint 


distribution of C, D and E, where C and D are non-negative, The transfor- 
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] be the joint moment-generating func- 
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APPENDIX IT 


THE ASYMPTOTIC EVALUATION OF A CLASS OF DETERMINANTS 


ARISING IN STATISTICS 


The following method of evaluation of a certain class of determin- 


ants was developed in detail McGregor [8 ]. 


In determining the joint moment-generating function of several 


quadratic forms in normal random variables, 


form is frequently encountered: 
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Example, For n= 11 and m= 2, the determinant has the form 
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Repeating this process (n-2m) times, we obtain 
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0 0 | 
—n,m =7,m=T a5 
where 
BT Par) 
: Pores (0) 
U = ° e ° 
=—m-T,m-T+1 
(0) d Pom=1 
B 
Then |Q*(1)|/= 1 (+ = 1,2,...,m-1). Successive premultiplication of | — 
y 
by Q*(1), Q*(2), ..., Q¥(m-1) yields 
B B 
_—_|= | Q*(m-1) Q¥(m-2) ... Q*(2) Q*(1) |_=_ 
F F 
B 
= tax ! (11.23) 


where F* is the mx 2m circulant type matrix with first row 
Sil 
(ike -N(1;1,3,...,2m-1), N(231,3,...,2m-1), eoog (bye N(m-131,3,...,;2m=1), 


(Ci) Nim; t,3,,. -4,2m=); 0,-omes.., 0}. 


By definition, 


N(431,3,.-.,2m-1) = i: a Sewer x Poy Poy ++ Pog 
1 2 


, j 
LS iy < i, Si b60.S i, <m 


Using relations (I1.8) this becomes 


1 
N(j31,3,...,2m-1) = ) . ° ° ° ° y D, DP, «2 -® 
915° +O, k k, k 


while 


al 


FES. EL) 
; = 7 i M , | Puiayyts 
a Yo Woy Jerkt dilw xivtamn aqya teblwouts 


wh ; LEY Ae See | ut (i-). i wae ar 


- Ol. - 


1 
NMS cr;cu-l) = —————— 
(m;1,3 ) 9195.0 
Let 
a, = (-1)4 ) ie ) DO. O.. 4 (II. 24) 
1<k_ <k " "2 " 
SS Bp oss 


Then Pi: Po» -++> @ are the roots of 


m m-1l m-2 
p +a, + a, t+... +a +a =0, (11.25) 
and 
, ais 
N(disl,g4etesem-l) = (-1)2 — : (11.26) 


m 


where we have taken a = 1. Furthermore, by equations (II-5).) G16) and 


(G82) 
om em-1 am-2 Tp 1 m 
Cee, Cow + Tac ae + cee, 
m=] 2 
eco ee +... + co? + 1? +1 
Dafpe i s _i 22 y eee 
= (9-9, )(9-9,) «-- (9-9, )(9 9, ey 9, 
ue m m-1 
a aa (Gp + cal ? DINE eel og sein a) 
m m=1 
x(a G +a 10 +eee + 8, 9+ 1) ee 
and, by identifying coefficients of o, we obtain 
Panes (a aera a ig jae as) Gii27,) 
j a, m-j Ll m-j+l ei jem 7 


(j = TSO tn) 


Thus F* becomes the mx 2m circulant type matrix with first row 


a a a a 
frjestel, patsé ot EO. 0} 
d a > a 3 LS Led) a > a B a > > > OF ONO) 


m m ™m. ™m m 


yd ,sronmetai0% 


tt 
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From equation (11.23), 


B 
a a a 1 
. a : (0) 
an an-1 *y 
B . 
peer = an (0) > 
EF 
he ee ig a 
and from equation (I1I.22), 
I i 
-In - 
K, =a (p,,_,7P ) | M(2t-1) | 
1 m i<i<j< m 2j-1 “2ai-1 i 
B 
a a a 1 
™m m-1 1 
(0) 
a an-l ay 1 
x 
(0) : : : ; . (11.28) 
a al: ay 1 
A similar reduction gives 
. 4 , 
Reo i. Apia cba,) | I (ae) | 
S  ldicjen 9 43 t=1 
af a a a 
1 m-1 m (0) 
1 a, eat an 
x » (11.29) 
(0) L di An-1 an 
H 


Using equations (I1.8) and (11.24) it may be verified that the 


following relations hold: 


fi 


| 
: 
| 


i 
= 
H 


em 
| 
ng 


ee eo eee 
5 eas | 
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en Ne (-1m(n-2m) y ~(n-2m) 


) ) ciate ) 
5: ol ly ; em mn ’ 


(-1)2n(m1) 


(p.-p,) 

I<i<jgom tJ 

ey (m-1) 

(-1 M+s5m fl ( Pp) 
Rei a 1 - 9,9.) Gueon)a 
oa Vier os ate 
m LS ed = 


em m m 

=] 2 yy 
I] M(t) = in 8 (1 - 9.9.) (ORT cei) sts 
tel aclencl) pills va lhe ¥ 


(Goes Boe =e), OO) 20.) 
ej-l 2i-l 1 i< jem j i 


w 4m(m-1) 
cy] (9, - 9) 


coe 1<i<j<m 


AD teal = 
i<igjen 7424. 


Finally, using these relations and equations (11.2), (11.15), (11.28) 


and (11,29), we obtain after simplification, 


m 
(0) -(n-2m) c -1 
D we el) (1 9.9.) 
i, j=l 
a ame MOAN Sol) ctl 
fm Ome a4 : fom Amd nae, 
(oO) | | (0) 
a a a 1 a a a il 
x m m-1 1 m m=1 1 
e ° re) 
(0) . : : (0) 
am of: CA a ; an mel i 


(II. 30) 


where J, is the kx k matrix with units in the secondary diagonal and 


zeros elsewhere, 


iieupe bas, prety 


i 
ji 


rot pashy. 
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APPENDIX IIT 
THE ASYMPTOTIC EVALUATION OF TWO DETERMINANTS 


ENCOUNTERED IN THIS THESIS 


A.) Known Means 
In order to obtain the joint moment-generating function of CG, D 
and E (2.6) the determinant of A has to be evaluated. A is the non- 


singular, 2n x 2n partitioned matrix, 


| 3 
A = |%1 1 35 (111.1) 


with the n xn_ submatrices 


1-2T “Py @) SS reatee) 0 
2 

“0, 1+p)-2T “Py © 0 

0 12 Lip, -20... 0 0 
Q, = : ; (111.2) 

O @) O So) oleae eae a) 

[ © ) @) =P, eee | 
1-28 Sons 0) SO) O 
=O 140--28 =O O O 
2 2 Des 
2 

O “Po 140,728... O @) 
Q, es ; (III. 3) 

0) © 0 ae rom= 25 “P.5 

fo) 0 0 “0, 1-28 | 


and Q, being a scalar matrix with diagonal elements (-U). Since the 


matrices are nonsingular, we may write (III.1) as 


eer bant 
i NAMEN ees % 


Qa ! Q5 &\ 2 9 | & 
-- + --- “ =) wi iy 
- 1 
iad Fi : = — - ; ‘ 
| | | aL | 
Cp iby CE Np yiche tee i ee oe 
Vou 
Bitlnd@| (ln, Lew lQay Qe 
2) Fs 2 5)ss84) ee esa (III.)) 
Thus from (III.1) and (III.4) the determinant of A is 
| eee 
my 3 (8 Bang ae eo 
CfA |e “feat ut (etl | See [eee ie ceed ae 
Ore os) ear ay ea 
| -1 
'Q,° Q 
= aoe (111.5) 


Again employing the theory of partitioned matrices, Aitken [1], consider 


where 


| 
FS 8 le 2 
Se eneeg Rie et wr ede ay viet oleate clara 
| 
| -1 -1 
= -1 Pe eer ay ee cs ee 
Siero ie 
= 3] Taree fo ce CRTC AEST 
IQ, Qs 9 1 Qs 
SF Es ae, | (111.6) 
— 1 5) 2) b) > ° 


8 is a scalar matrix. 


A, to saselereted 


i 
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Substituting (III.6) in (III.5), we get 


-l ~-l .2 
|A| \. IQ 1IQ511z a a) Q., Qs | 
2 
= |9,9,-Q5l - (111.7) 


Thus using the notation of Appendix II, we may write |A| as 


1 XX Xy LO ee ONO 9) ) 0) 
1 x Cee Meee: Om LO 0 0 
|Al=p465 eer (EEES) 
(oj) {0) 0 Oo OM Oe eck x) Xp x1 1 
0 0 0 0 0 0 e688 0 Bay Bis Bio Bay 
L- 1-25 y2 
where p., = +1-—— ‘ 
11 Py Pr P1P5 
2 
_ =(1-2T (1495-28 ) 
pw alice) 
1+0, -2T) 
-(1-25) zn ( 1 
hl iS p , 
2 1 
(1402-2) (1492-28) ais, 
-(1+9,°2T 149,28 
% s 1 = Basi ’ 


(1407 -21)( 1405 ~28) y2 
a + a Sees 


Bis = Boy, = l, Bay, = O, Boo = X, and Bos # x): 


(Note; with regard to the notation used in Appendix II, me= 2 in (III,8)). 


Following the method of 


4} x } determinant, 


where 


nu 
[Al = e,°, FEET WC 
Seen ee tetera) ta 
Boy Boo Boz Pon 
Bon Box Bop Bay 
H = ’ 
fe Pi Pas) Pye ef a7 
n-l n-1 n-1 n-l 
ae re ee re 
P) (n-4) = ; : ? : ’ 
SH: Ba oe 
n-) n-k n-) n-) 
Py Pp Pz Py 
AU Neon utc ta 
Bore bs Be 
a Cee fea es Ph 
Bh (O)= (x, .)= > 
By pea aes. Pi 
rides ae gai me 


and Py» Po» P 


By) 


oe 
Pp 1P 


3 


+ X 


BP, (n-h) By (0) 


2 
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Appendix II, we may ultimately reduce 


; ah are the roots of 


p- + X)P + 1=0 


(22U.11) 


From the symmetry of the coefficients in (III.12), it is clear that we may 


write 


a 


= 


| ) 
: 

mf 
be 
| 
b 


am oy Jott) 2 


- 1Ol ~ 


> 1, = + 0 = and a4 (ILL.13. 
vy} = Py Po ees —,” Ps ~ ? © Py cs} 7) > N Waar yp 
J 2 
and by [11] 
QM, + L + QP. + Be -xX 


(1404-27) (1405-28) 
OSH abisak ot As 
2 


(TTT 1h) 
and (IIL. 14) 


1 
tal 
i 
nv) 


1 1 
(9, + 9, ) (2 + 95) = X, 
(1495 -2T)( 1405-28) 
P40, 24°, 


ran 


The transformation (III.14) may be considered as consisting of the success- 
ive transformations 


(1404-21) (1405-28) 


2(P+Q) = Ls ; 
eal G2 
(itd, 15)) 
2 2 
aes (1+) -2T) (1405-28 ) i 5 3 \F 
S GuPRei aon CMR anon, iw Tae 
ir 2 12 il 2 
v(r,+r,,)-Tr,,-Sr i 
(where we have put U = 2a eal », and ve=O0 following Appendix I) 
Tyo 
and 
i 
Pliee  h=meb 4 
9) ®, 
: (III. 16) 
tie ec 
Po ®, 


Transformation (III.15) will map paths in the T and S_ planes, parallel 
to the imaginary axes, into paths in the P and Q planes which,by the 
proper choice of branches in transformation (II1I.16), will be mapped into 
the interior of the unit circles of the ?1 and 0 planes, Thus we may 


take 19,1 <ealleeneniicl lo,| < 1 except in the critical cases where equality 


; SOW wh lagi SIE A UD . 
abelian: tnan emshreteetio ners nalntiae aie cbdosnenitacead 
ft “D vie 


ae 


}: 
} 
| 
y 
id 
ma 
i 
Mf 
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will hold, 


Therefore, since each element of the original determinant |A| is of bound- 
ed modulus as n-+o (see (III.8) and (LII.9)), |?, | <li: IP, | IS 


and «, are independent of n (see (III.9)), and m is 0O(1) with respect 


to n, we are justified in using the method of Appendix II to approximate 


(0) -1 
A , re ayn o - ° 
typical element of B rn Bi (n ) Py. CO) any (EID. 10); as 


Since p, = ®1> 19, | < and) sp 


° 5 n | 
where the terms omitted are relatively O(e€ ) for some |e; < 1. Hence we 


may write 


BB, (n-h) By (0) 


B n-l n-1 
ve Po Pi Mo Neo "os Tol 
hue is | bo (GhED 7.) 
Po Py ed me ho” Py s hip 
n- n- 
Po Py, 
n-) n- 
Ah 


Thus, using (11.15) and substituting (111.17) in (111.10), |A| may be 


approximated by 


4 


| ‘Bow 


i 


oh ‘peeban)'| 


oe 
Gn 


re 


fa 


Ra eon aa) a) 


X11 Mo M13 My ly Myo yz My P (111.18) 


JOT ire A 


The elements, x of By (0) are of the following form (see Aitken [1 ]) 


alg)" 
(n1y> > Ty) 3? m1) 


1 
 (Py=Pg)(Py-P5)(P Py [1,-(Po+Pz+P)) »PoPztPoP),+P5P),2~PoPsP),] » 


(s054 Moo >M oz 2Toy,) 
= [1,-(p,+P,+P),)»P.PotP,P) +P P) »-P, PsP) | 
(pepe a PnP) Le ake ak 
and similarly for COLE ur turie (my, >My o9My 59 Xp) 5 


Then, following the reduction procedure of Appendix II, the two determinants 


in (III.18) become 


1 
~ (Py-Pg)(Py-P3)(P~Py)(P5-Po)(P5-Py,) 


B 


. 1 -(Po+Py) PoP), 0 
0 1 -(PotP),) PoP), 


and 


ig 


Halil | 


Ci me P . 
ae rac 


‘ 


ie 
= 
=~ 
. : 
=e 7 " Te 
« (Mae é 


add tL atbneqaa do: 


Se Pa WAS ey F 
Nh NWSE rapa Newey kre btm are 


Hine 7 i, f 
a mT 
i ey ; ey 

4 a9) 


i 
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1 “Sane neneeieeeaiamneneemiene ean 
Ty Myo Myx Nyy (P-P,)(P5-P3) (Py -P,)(P,,-P 2) (P),-P 2) 


H 


x| 0 1 “(pi+P,)  PyPs 


H 


Using the transformation (11.2), we get 
a, -(9,+9,) and a, = 9) % 3 (III.19) 


and, from (III.13), we have 


© 


Soy Ges (pe) ee eM pp. =e and <(p.+p,) = 2, 
ak a, ? ey, a ; 1°3 2 Maas aes 


Thus (I1I.18) becomes 


|A| ~ 0105 { (-1)9;0""*) [(P,-Pp)(Pp-P5)(P5-P,)(P,-P,)] 


; ge dul a ff ay m5 O 
ene (IIL.20° 
as a, a, i! 0 Orn d a, 4, fe (IIL. 20) 
fe) a, ay 1 H 


t 


Expressing (P1-P5)(P5-Ps)(P5-Py,) (Py -Py,) in terms of ay and ane we have 


L iv 
(p)-P2)(P5-P),) =(9, - 9, 25) 
1 2 1 
= 9195 - aay (P +0,)° + 2+ 
per RR oe 1? 5 
1 2 1 
= ances lan) a2 ae 
2 ay 1 a, 


1 
a pes) asia 


Ve ented: 


and 


| 
i 
=2- 9,9 
LM SO Ms 
1L 
= 2 = re . 
2 a, 
elk 2 
=> Gt a, 
Thus 
oO" oO 
ee) 
Ale n-4 2 
a, (1-a,+a,,)(1+a,+a,)(1-a,) 
En figa!° eet ae 
Boy RX 1 O IL ay 
x 
a, a, i O 1 XX 
O a, ay uy O i Bin 
Seem gy aah eae eye OP Boa a sao ene Ue 
From (111.14) and (111.19) 
91+ 
ie) 
x= - (9, + % + ] 
1 1 2 ? P5 
a 
1 
Sj (1 +a ) 
a, 2 
and 
2 2 
Kags DoD ie mn + 2 
Be DU or wm PND 
di 2 2 
= a, + 2, + a5) 


Finally, using 


(III.22) in (I1II.21), we have 


3 


(III.21) 


Hel norstanne 


g (, 


ilidt 


7 te roa. NAY 
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n 2 : 2 
Py [a5(By4785) - 85(By 578, ) (B51 744) J 


2 . ; 
(1-a, +a,,)(1+a,+a,,)(1-a,) 


where terms which become exponentially small have been omitted. 
B.) Fitted Means 


In order to obtain the joint moment-generating function of C*, D* 
and E* ( 3.7) the determinant of B must be evaluated, The approximate value 
of |B| is the key to the problem of finding the distribution of the projuct- 


moment correlation for the fitted means case, 


B is the nonsingular, @nx2n, partitioned matrix, 


R R 
ae (III.2h) 
with nxn _ submatrices 
at+a “P, +e a eoe (oF a 
“04a b+a, “P40 S60 a a 
a “Py +a b+a, esis a a 
Ri = e e ° ° e ° e e ° ° e ° ° e ° e ° ° e e ° 0 ° > ( Til ° 25 ) 
a a Qa ooo b+a -0 pie 
a Qa Qa eee ad @) 1 a+aq 
al 
"| 
c+B -P+B B Jb B B 
“Po+B d+ -P.5+8 ° B B 
B -PAtB d+B cee B B 
LIL. 26 
E, °° e J e °° e ° e e ° e e e e e e ° e e ( ) 
B rf) B See d+ -“Po+B 


a, kee mood t 


i 


= 
i eye ay: 
Hes eo 
ahd bY Des iT ee ery tei 
aU au 38 7 RD, De. sae oats puthaat 40 more 


$e ro) 


ae 


=) RETA MAY Dono k Abas we | ee 


gi! ah 
ae ee ra 
+ f 
Fe $ 
ii 


ve ey AMIE SN PTT 


oe Tees 


and 


where 


and 
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7 Sed 7 Y 
Y care y Of 

-U+7 oe i/ Y 
aM R ae GaaMel ee te ; (ELE 27) 
7 or -U+7 Y 
Y sieite 7, ~U+7 


(111.28) 


In order to be able to use the determinant-evaluation method of 


Appendix II we must transform |B| to the form (II.1). To do this we shall 


use the bordering technique of Aitken [1] and Muir [10] followed by a 


series of pre- and post-multiplications, 


The determinant 


determinant, 


|B| may be written as the (2n+2)x(2n+2) bordered 


ms Cert eee == 


wa 


ootteu love: sasmtmgesab. ef 


or. ue 3) arto! 


» yo bewollot [or] thuM’ bas [ 


1 ha were aca laud Tinie 
0 0 
ln | 45 
Ic l= a9 0 (III.29) 
0 0 
- | OR R 
a5 =2 
0 0 
Oh le Yow Ta ae ose eh | Cli, a OLY. aiioteveg: TCG | 1 


Premultiplication of GC, by the (2n+2)x(2n+2) partitioned matrix 


1 
Us Or a0 eee OF 40120) £0 eee (OA OF Rao) 
-1 1 ce) eee @) 0 eee (6) 1) 
-1 @) i eee 9) O (0) 0) eee (9) 0) @) 


yields, 


em benokatttaq: (Saad) xf 


+08 ii i ats. x 


anit richie et renee eam ants aa arama ts ARMA ri damental eda rt 


Oa 


! 09 = 


Lelia a a : a a vi 7 Y 
-l a Py O O @) O -U O 
-1 “P) b “Py 0) (9) 6) 0) -U 
-l @) “Py b 6 (0) @) 0) O O 
-1 0) O 0) a b “Py 0) O O 
-l (@) (@) O coe "Py b “Py O O 
-1 i O @) @) ° ° 0) a) a 6) 6) 

(@) -U (@) (@) e oo 6) 6) 6) Cc ao 

@) (6) -U 9) ahone @) O (0) “Pp d 

(6) (0) (@) -U owe ) (6) (9) 0 “Oy 

; | ° e ° e ° ° ° @ ° e o ° ° ° ° ° o 

(@) (9) (@) (9) ere =U) (9) O 0) 6) 

O 0) 0) 6) eoo° (6) -U O 6) 6) 

(9) (9) (@) O silexe (6) (@) -U 6) 6) 

) Vf of 7, Somme 2 Y B B 


which we may write as 


rt | Qi cere nO | y AG 
-1 
; Q a; 
Io. E 
C = as aT 
—2 0 
5 Q 
0 | 3 42 
Oe cco 7 | ee 


~ 


BD 


| O---0 jo 


i 
re 


i] 
Rico e 


= | 


° if Tey aes Aoi 
O O 0) 0 
° O O 0) te) 
O O 0) ) 


a d “Pp 0 j-l 
Rie “2,5 d “Pp -1 
° (@) Te) c -1 


(III. 30) 


5 (III. 31) 


where Q, and Q, are defined by (III,2) and (111.3), respectively, 


“Ha 0 ae. 


Fler opicn Hat Rorelueun sin vet seas Raubisnt yada Alene dst¥! ell Dei Rinse nesab ernst heen ea tak Rice en Dede sy 9h 


G: gf” 
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and Qs is a scaler (nxn) matrix with diagonal elements (-U). 


Premultiplying Cc, by the (2n+2)X{2n42) matrix 


i 
H 


we get 
Dae ceara teins eae ay eee wai Fac 
=] @) 
21 2 
le. j=|z.c,| < e: | ® (111.32) 
=% =2=02 — SSS SSS SSS = , 4 
0 =i 
i gs a, | 3 
0 =1 
oe ee ee ee 
£ |KO gud. A, Os\Behe tom. seme cen 
and premultiplying c. by the (2n+2)x(00;2) matrix 
nwt a = 1 
1 Ou gas 0.) (Bo ] 
Te tate at ger 0 
1, = 
iB) °. eo ° ° ° oo o °o o ° ° e ° eo o 
O 0 ooo 1 0) 
oO 6) eo°o oO 1 


we have 


Sheen er atietgw erin morn r meyers hareeg 


ee 


pied ons ab: i 
crotondg Ba | 
ay 8h) 0-8 


‘= 
ne 


Xs | « ens (710 | Ol ate (20a xy 

oil @) 

al a, 

4 0 iY 

Ic, | = |2,¢.| = aay eee ema oR eT aT ce ’ (ILI. 33) 

3, 2 

: i 

crue gona, | Oi | ee ae 

where 
dee eas 2B -E , 


5) 4 ap-7 ; fat 


To procure more zero elements in the first and last rows of Ic, |. we proceed 


as follows: 


firstly, solving for A. and A. in the two equations 


i 2 
w~ A, (b-20, ) - K3(-0) -=§ O 
and 
Be A (-O) ae A,(d-20,,) = 8) P) 
we get 
Nat, i ofS Poe) 
3 V(d- Ye 
[(b-2p, )(d-20,,) -U=] 
and 


Xr = . ol a 
2 [(b-20, )(d-20,,) -U=] 


aarp mi ake albany si ents Ate. dinette 


EG ake ae ee 


| ; at 
esl Boe teekt ed? ab 


Sy 


‘ot 


‘enotisepe, ows 


Then premultipl 


1 -A 

O 1 

O O 
iy 

1) 0 

0) O 

0) 1) 
gives 

C 

Ics | 


lle - 


ication of cy , (111.33) by the (2n+2)x(2n+2) matrix 
1 sty sh Ne re “5 “Ay Y 
Oe @) OF mci oe. 0. 6) O 
1 5 e) O O 0 16) 
(Os G0) epee eae) 0) 
(9) °° 0 0) oe ° oO 1 (6) 
(@) oo ) 0 eae 16) 6) 1 
1 
ms “2 ih By 2 
= |Z,.¢, | = [(b 20, )(d 20) U=] 
Xe | Xe (sisi en) Xe | a OmmeaecucO Xo | Xe 
= 0 
: Qi Q5 : 
il 6) 
| = ) 
0 | =i 
0 | gs; 5 “1 
xy | 0) (OMe ey ua 6) 0 X, x a x5 Xp | 1 
(LIL. 34) 


where [(b-20, )(d-20,.) -y2y"! is factored from the first row and where 


be 
iH 


[(b-2p, )(d-20,,) “UII x, + nay] 


x,{(b-20, )(d-205) -U?] + no,(d-20,,) ; 


a ity 


Ate iy 
. 7 Pe yo a 


= iis 
Xe = [(b-2p, )(d-20,,)-U][a-A, (a0, ) + Av] 
. af (b-2p, )(d-20,,)-U*] - a(d-2p,,)(a-p,) + a= 
= a(d-20,,)[b~20, -a+0,] 


= ap, (0 ,-1)(d-20,,) 5 


% = [(b-20,)(d-20,)-US][A,U = Ap (c-,)] 
2 a(d-29,)U - a(e~p,) 
=o U[d-20,-c+0,,] 
= oP .(P,°1)U 


and 


Xe = [(b-20, )(d=20,,)-U*) [x,, + nn, 
= x, ((b-20, )(d~20,)-U#] ¢navU . 
Secondly, from the two equations 
Xo Ng (4-205) = A, (-U) = 0 
and 
= Ag(-0) “ A, (6-20, ) 20 
we have that 
nN %,(b=26 ) 
3” [(b-20,) (4-20) -U4] 


and 
x,.U 


N° Tbe eee, 
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Hence, premultiplication of C_ (III.34) by the (2n42)(204+2) matrix 


5 
ul 0) 6) fe) fe) . O 
0 L 0) ° 0 0) co 0) 
0) O 1 ° 6) 0 ° 0 


gives oc. ; 


where 


(@) 0 O one (6) O Ayays (@) 
fe) “Ay “hy Bo “hy, “A, aig A, 
“2 
ISe| £58. | Lo [(b-20, )(d-20,,)-U*] 
Xe | Xe (Vane pease (a) Xe | Ol vevie 
«1 
° Q, a 
ih 
n corm eee: | ores ate towne VE BON A arene eee ereienere ese: 

(9) 

i Qs S5 
Xo | Hig 0 cee O Ry] Ry 0 eee 


cite | eae 
0 
0 

oe ee 
=i 

et al) 
(L111. 35) 


[(b-20, )(d-20,,)-02]"? is factored from the last row, and where 


Xq = [(b-20))(d-20,)-U"][nr, + x, ] 


= (are es, We 
= nxU + x, [Cb 2p, )(d 20) U ] 9 


Pra ' () 
sa ¢ xe 7 a td 


Rha laieahiaeywinn | ke 9 (oer en CARMEN NER g AUR Ihd nia pf ni Na © sad ated HCH pe ay ReaD WAR da pm afi ildinalna is 


y (} 


et Ra (a nO i tena le arama er hd EL Se Ce ANN p= Paine 


Bt banc MeL IY  ARe RUE ER pn en ya NAB Gap 


Z 3 és ¥ Oy : in ert 
i he ps 


Woy steal a 
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ie [(b-20, )(d-29,.)-U2)[~A, (a~, )4A,0] 


= -x,U(a-p,) + x,,(b-2p , )U 
= x,U[b-20, -a + 0 


re) 


= xP, (0,-1)U ’ 


x11 = [(b-20, )(d-20,)-U"][A,U-A,(c-,,)+x,] 


= x,U°-x,,(b-20, )(c-0,,)+x,,[(b-20, )(d-20,,)-U?] 


- x,(b-2p,)[d-20,-c+0,,] 


= X,0,(0,-1)(b-2p, ) 
and 


_ [(b-2p) )(d-20,)-U7][nA,+1] 


= nx,(b-20, )+[(b-20, )(d-20,,)-07] ‘ 


To replace most of the (-1)'s by zeros in the first and last columns of 
ic,|, we shall employ the same technique as in the derivation of Iecl, 


(111.35), from Ic, |, (£11.33). Solving the equations 


«1 


i 
GS 


H, (b-20,) - un(-0) 


and 


i 
So 


ae (en ay) zs 
Hy {-U) - u,(4-20,) 


we get 
-(d-26,,) 
Sey [(b-20, )(d-20,,)-U*] 


and 


| ae i ll ai 


i" j } \ i i { 4, i 
’ | ; ) | | | iu ne y bad ¥ Biche ae { 
Dement | | | (bapfaal (thud gai-b){ 
y ‘ = F i nd . k Min | ay i _ Ra aay 
i rd 
See RI i AP tT ole Hy 
> to enmutos Jeet bde Jeukh 449 at actor ve 


SOE AD a RUE 
saat iiisdidas 


Aide eups 


hE 


Ait va) 


M2 ~ [(b-20,)(d-20,)-U) * 


2 


Post-multiplying C- (111.35) by the (2n+2)x(2n+2) matrix 


=6 


the determinant of the product becomes 


IS| = IG.t,| = {(b-20,)(4-20,) -u2}? 
= O AOS ) x x O stone O 
Q Q, 
Q, Q 
| 


(IIL. 36) 


keen | CSqinst 


ji a) 2 : My j 
ded Jaubo: 


ve 


‘gemtoned Jo 


where [(b-20, )(d-20,)-v2]"" is factored from the first column and so 


Y= [(b-20, )(d-20,,)-U"][x,-1, (2x,)-u,(2x,)] 


il 


x, [(b~20, )(d-20,,)-U*] n 2[(d-20,,) x, + Ux, ] 5 


va = [(b-2p, )(d-20,,)-U"][-1-n, (a-0, ) “5 HU] 


= -[(b-20, )(d-20,,)-U®] + (d-20,)(a-0, )-U* 


>) 


es -(d-29,,) [b-20 -at+p, ] 


1 


= ~0, (0-1) (d-20,) ; 


= [(b-20, )(d-20,)-U*][u,U-u,(c-0,)] 
= -(d-2p,,)U + (c-0,)U 


= -0,(0,-1)U 
and 


-Yy = [(b-2p))(4-20,)-U"]Exg-H, (2x15) - Hp (2x, 4) 


] 


Xo[(b-20, )(d-20,)-U] + 2[(d-20,)x,, + Ux,, 


Similarly for the last column, solving the equations 
“uz(-U) - u,(b-20,) = 0 
-1 “uz(d-205) = y),(-U) =O ,; 


we obtain 


-(b-20, ) 
M5 ~ [(b-20,)(d-20,,)-07] 


 amokteups edly 3 


yeed! 1, 
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and 
u - U 
 ~ [(b-20, )(d-20,,)-U=] 
Hence, postmultiplying c, (111.33) by the (2n+2)x (2n 
1 0) fe) So 9) fe) ce) 
(@) 1 (9) eae 0 (@) “Hy 
0) 1@) 1 eee 0 0 eat 


_|F 
oe) 
oO 
e) 
Oo 
(o) 
q 
r= 


Oo 
Oo 
Oo 
(oe) 
fe 


io} 
oO 
oO 
‘e) 
(@) 
im 


we have, as the determinant of the product, 


CL| = [(b-20,)(4-20,)-02)"" 


lee OO cas 4 0 4 & |x, Ore, LO 


to 
— 
we 


2) matrix 


io 


oh ll 
Wye 


~renn cm nonersoaraclraiini tit Noy 0 10h bay repens meh eh oh 
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where [(b-20, )(d-20,,)-v2}"* is factored from the last column and where 


5 = [(b-20, )(d-20,,)-U"}[-H, (2x, )-H,( 2x, ) + XQ] 
- 2[Uxc+(b-20, )x_ ] + x,[(b-20, )(d-2p,,)-U=] 
¥¢ = [(b-20,)(d-20,)-U*][-u,(a-0,) + p50] 
= U(a-o, )-(b-2p, )U 
= -0,(0,-1)U ; 
Yq = [(b-2p, )(d-20,,)-U*] [uy “Hz(e-e,)-1] 


-U> + (b-20, )(c-,,)-[(b-20, )(d-20,,)-U*] 


-(b-20, )[d-20,,-c+0,,] 


-0.(0,-1)(b-20, ) 
and 


| Yg = [(b-20, )(d-20,-U"}E-1y (2x, 4)-H,(2x, 5) + Xo] 


2[Ux, )+(b-20, )x, ,] + Xj ol (b-20, )(d-20,,)-U*] 


prays 
pannus 


He 


i 


} i i 


where, using 


11 


relations (III,28), 


= %,0,(0,-1)(b-20,) = (8 - 


In| = [Gp 


ap, (e,-1)(d-20,) 
a,,(P,,-1)U 


xP, (0, -1)U = (Bp 


it 
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[(b-20, )(d-20,,)-U2]" 


6 | 0 fe) 
| 
ion ° 
ale) 
i 
XE); 
Te i 
for )e4 (0, 10 
1 
Oz) > 


o(4) 


2 
= )e,(9,-1)(b-20, ) 


Ve ii} 
; i) L 
{pie ARV iNew ele a magne lata nae mali 


tenia tte ehartelinnepnaeh aban eatinetrieae  T 


yy 


46 


ti 


x, [(b-20, )(d-20,)-V*] + 2[(d-20, )xc+Ux, | 
x5{(b-20,)(d-20,)-U2]* + na(d-29,,){(b-29, )(d-20,,)-U] 
+ o(+) 
laa [(b-20, )(d-20,,)-U2}° 
+ 2T(d-20,,)[(b-20, )(d-20,)-U*] + 0(+) ; 
-0,(0,-1)(d-20,) , 
-P,(P,-1)U , 
Xo[ (b-20, )(d-20,)-U*] + 2[(d-20,)x,, a Ux, ,] 
nx,U[(b-20, )(d-20,,)-U*] a x, [(b-20,)(d-20,.)-u2]* 
+ o(2) 
een | U[(b-2p, )(d-20,,)-U*]- oT [(b-29, )(4-20,,)-U2]° 
+ 0(2) , 
2[Ux,+(b-20, )x_] + x,[(b-29, )(d-20,,)-U*] 
o(4) + x, [(b-20 , )(d-20,,)-U2]° 
+ naU[ (b-20, )(d-2p,,)-U*] 
Ea [(b-20, )(4-20,,)-U2]* 
+ aTU[(b-20, )(d-20,,)-U*] 4s o(=) j 
-0,(9,-1)U : 


-P,(0,-1)(b-20, ) 


(III. 49) 


(III. 40) 


(III, 41) 


Com Come 


and 


il 


Yg 2[Ux, ,+(b-20, )x,,] n X1ol (b-20, )(d-20,,)-U*] 


il 


0(=) + nx,(b-20,)[(b-20, )(a-20,)-U2] 


+ [(b-20,)(d-20,.)-U2]* 


led (b-20, ){ (b-20, )(d-20,)-U2I 


+ [(b-29, ) (d-20,,)-U2]° + 0(=) (112.42) 


We now expand the (2n+2) determinant of (111.38), which we denote 
by 


By = [(b-20,)(d-20,)-u2}* |B]. (112.43) 


Ignoring terms which are o(=) ’ 


B,~y, 3, + Cie Ye BS 
that is 
BLY yy a Ye Bs 5 (LIT. 4) 
where 
%6 
) 
Q 8; 
) 
%6 
cn 
fe) 
Q, Q, 
O 
ua 
10 Oo Goa © X16 1 fe) Bre) 11 Ye 


aN 
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and 


75 


yy X10 0 eee (@) x x 


Expanding B, and B.. and ignoring terms which are o(4), we obtain 


2 
Bo~ Yq IAl (III.45) 
and 
Bs mY), (Alas (III. 46) 
where |All , 
ore: | 8; 
a | % 


is the same determinant as encountered in the "Known Means" case (III,1). 


Thus we have from (LII.43), (1II.44), (111.45) and (111.46) that 
-4 
é _ 2 
|B] ~ [(b-20,)(d-20,,)-U*] 


«Ly 1¥e7¥5¥),] LAr (111.47) 


t 
ore 


oie Hobe ametad ry 
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Using the relations (II1I.39), (1II.4O), (I1II.41) and (III.42), we may writ 
4TS 212 
Yg-Y5y, = | [rast>| [(b-2p, )(d-20,.)-U=] 
+ 2n( 4-20, ){(b-20,)(d-20,)-U2)} 
eo (b-20, )[ (b-20, )(d-20,,)-U*] 


+ [(b-20,)(a-20,)-021°| 


- lise [(b-20, )(d-20,,)-U2]° 


+ etl (b-20,)(d-20,)-v4)} 


{fas U[(b-20, )(d-20,,)-U*] 
: + {(b-20,)(a-20,)-021°| + 0(5) 
-u2]* 


= [(b-2p, )(d-20,) 


28(b-20, )[(b-20, )(d-2,,)-U?] 
Ws Be 
+ Es | [(b-20, )(d-20,,)-U ] 


+ [4TS-U*] (d-2p,,)(b-20, ) 


+ 


aT (d-2p,,)[(b-2, )(d-20,,)-U*] 
+ U?[(b-20, )(d-20,,)-U*] 

US d 2)- 
- Fed [(b-20, )( -20,,)-U 
- [urs-u2]u- 


+ v?t(b-20,)(4-20,)-U2]} sf 0(4) 


continued 


oh Ay 


S1u-(agvb)(jasea)] 


Palse of a 
E U~( a9 
ba 


dah 


-b)( ada) Ju | 


“pies | — SSt=( 08-8) 08-4): 


ree ie Fa t thy 
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= [(b-20, )(d- 29,)-U?]? 
‘jee [(b-20, )(d-20,,)-U*] + 28(b-2p, ) 
+ at (d-2p,,) + 20° -| Fave |l(>-20,)(2-20,.)-07) 


+ (uns-v2)} + 0(4) 


[(b-29, )(4-20,,)-U]? 
«{(b-20, )[4-20,5428 + at d-20,428]} + o(4) 


[ (b-20, )(d-20,,) U2)? 


i} 


1 
*(d-20,+28 ][b-20,+2T] + 0() 


and using the relations (III,28) 


Yig-V5Vy [(b-20, )(4-20,,)-U7]? 


2 2 i 
x [1+05-28-20,+28 ][ 1+0}-2T-20,+2T] + 0(=) 


i 


[(b-20, )(d-20,,)-U21°(1-0,,)°(1-0,)° + o(2) 
Thus (111.47) becomes 
|B] ~ [(b-20,)(d-20,)-V2]"* (1-0,)*(1-05) [Al - (111.48) 


From (III.9) and (III.22) we have the relations 


( 1405-22) ( 1405-28) v2 


fe) 


1 2 2 
= (l+at+ar) = 24 


and (111.49) 
2 2 
a (1+0)-2T) (1405-28) 


) 
2 Pa Po 


eRe ‘ re ea aa 
(“)O + i[ 884 nS+bjts + [ass 


: 
by * 


re 


i as~ag 


11884, 


peas raitey 3 Oe 
Jo + [1+ agets~ Foe 
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Now 


[(b-2, )(d-20,,)-U*] 


= bd + 20,0,-U® - 2(p,4+0,.b) + 2p 


5 io 2 


so that using (III.28), 


[(b-20, )(d-20,,)-U*] 


(1405-2) (1405-28) 


Uz 
= —P.A —_—oO HH Oe 


eee P45 


(1405-28) (1404-27) 
-2 |——2—_— 4 +2 


Pn Py 
and substituting the relations (III.49), we obtain 


[(b-20, )(d-20,,)-U*] 


2a 
= 0,0 las (Uyeceas) op (l¢e.) + 4 
ay ay 2 


l2 Dae 
9,0 
a2 [ eure. 
= a, leaytast2a,+2a,a,+2a, | 
P.0 
ae (l+a,+a stan 
a, Lee 


Finally, substituting (111.23) and (I1I.50) in (II1I.48), we have 


(1-0, )°(1-0,)* a, 


~ 


2 
0Po(1+a,+4,) 


—— 


nn 2 2 
_ PiPraln(By1745)-95(8y 9744) (854741) 


n 2 
a, (1-a,+a,,)(1+a,+a,)(1-a,) 


or 


ree, g . 


#f)f e+" us + 


U 


anode ov ,(Q4. 177) snolasies od 
te , D MoE a ay, 


re a | + eye Ty ea 


et es aK 
BS +, BS+ AS 


a real Sy 2 ( eet ca 


n«l n=l (yc 2 
P, 5 (1 p,) (1-05) 
= n-l 

a 


2 


2 i : 2 
: [a5(B)1-49)-45(8 15741 )(Bo4-4,) (111,51) 
ar regu ECE , 11.51) 
(1-a,+a,)(1+a,+a,) (1-a,) 


where Buy? Bro and Bo, are given by (III.9), a, and a, by (TET,hO) 


and the error is 0(=) : 


Kc) 


ee mah ie dean went ol ein 


a 
eh. 


APPENDIX IV 
EVALUATION OF AN INTEGRAL ENCOUNTERED IN THE 


RENORMALISATION OF THE DENSITY FUNCTION OF r 


Consider the integral J(m,k,a), 


2\m 


2 aod lf 
l-r lea -hor 2 + Lt, ma 
xr b 


J(m,k,a) = 
JO (£(tea)? hor? J241-01) (140)? boar J? 


wheres -lGn<ait, qdajecl, m< k 4 


Let 
(eaale ler Ai(tesag lee | (1v.1) 
Then 
lar = (ieee - Gaao)- 
= [l+a-(1+a9) ] [1+0+(1+09) ] 
= a(1-)[2+a(1+9) J 
= 2a(1-9)[1+ S(1+9)] , 
Poe S(1-9)(1 + ni 1+9) ] 
and 
1 1 a 
r = fa(l-9)? elves i l+9) 1 5 (Iv.2) 
Also if 
(ee) sana] (Giea)a 


l+2,+a - tor = 1+ 29+ ation 
Be) ep dinee = alice + 29 - al 


high = Nee = aila@ een oe) 


AA Lee i 
a naval Ni 


at Wk aAasTMUD 


hi ey ‘A 
(0 


Wa(L-r@) = af2(1ep) - o(1+p)(1-0)] 
= 2Qa(l+p)[1 - (1-9) ] ’ 
then 
ter = 5 (1) (1 - $ (1-9)] (1Vv.3) 
2 2,5 
(((l+a)* - far“J2 + 1 - a} = (leap) + l-a 
ela A o)) (IV. 4) 
2 2,44 
Gidsa)= see je 4 14 a) = (sep) + 1+ 
= 2[1 + ai 1+9)] (1V.5) 
and 
Ga) etter Ie = Cece) (1v.6) 
From (IV.1) and (IV.2), we have 
(L4a)* - har® = (14ap)® , 
-Sardr = 2(l+ap)a dp , 
dr = - ae. dep 
and so 
-(1+ 
dr = T T (IV.7) 
2 2(1-9)?[ 14 (1+) ]? 
Also at r =O, (1a)? = (1ecm)* so geal, 
and at r=l, (14a)* - lg = (leap )* , (Iv.8) 
(ea) = (eam) eigcim = =1 


Now substituting (IV.3), (IV.4), (1V.5), (1V.6), (IV.7) and (IV.8) in 


J(m,k,a) we get 
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1 4 
(149)"[1- S(1-q) 1” 22[14 Bi 1L+p) 1° (1400) 


1 
J(m,k,a) = hi aaa. a Ee de 
J 2™ aki 1- &(1-9)]* (14am) 2%°(1-9)2(14 ¥(149) 1? 


-1 
i m 
Sra “= k-m = (Iv. 9) 
2 JG IS Ue G.ca)e 
Let p= a -1, theny = 5(149), dp = ab 
l1-gm= 2(1v), Ll+om= oy 
and at g=-l, pe=O 
at Oban el, V = 1 
With this substitution (IV.9) becomes 
i m,m 
1 2 2 
J(mk,a) = “a7 k-m 43 3 
2 Toe etl ean) 
iL m 
== fi +, ; (Iv.10) 
eI (lene) aa) 
Let AG =1- WV > d€ 5 -dwy 
then at w-=0, €=1 
and at wv = ills C =O, 


Thus (IV.10) becomes 


1 1 
(Gee cae (ae Ge ™) at 


1 


By the theory of the hypergeometric function (see, Rainville [16] Thm. 16) 


[(ec 


: b-1 c-b-1] -a 
r'(b) T'(c-b) | Car cist) (1-tz) " dt , 


F(a,b;c3z) = 


if |z| <1 and Re(c) > Re(b) >0 


a 


Ee MM Senalepdonaiietiter sitet senate pa. 


Opt) 


ee) 


Thus 
il C(AYE (mel 
JI(m,k,a) = yes = ; 3 
ave ['(m + 4) 
2 
= a7 (k+3) 
Consider the case k = 2m + 2 4 
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F[k-m, = sm + 3 ec 


ees) ral semstese m + 


2 3 3 at) 


Then (IV.11) is of form 


JI(m,2m + : » a) = a7 (2m+2) B(m+1,5) F(m + A 5 ; m+ 2 

By a well known result [Rainville [16], Thm. 21] 

F(a,b;c,z) = Gest F(c-a,c-b3;c3z) , 
we see that 

F(m + 2 = m + z a) = (eeyre F(O, m+ m + a a) 

& Genre ; 

Thus 

_JI(m,2m + Ae a) = 27 (2me2) 15 Gaiety 1) (1-0)? ; 


(IV.11) 


(Gaevle) 


_ (elt 
oliivat 


J ond i [a ig) 


uy 


ray Cord Ve ite if oe ot} fl 


Wt 
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